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Abstract. We present a simple hybrid dynamical model as a tool to investi- 
gate behavioral strategies based on trend following. The multiplicative sym- 
bolic dynamics are generated using a lognormal diffusion model for the at-the- 
money implied volatility term structure. Thus, are model exploits information 
from derivative markets to obtain qualititative properties of the return dis- 
tribution for the underlier. We apply our model to the JPY-USD exchange 
rate and the corresponding lmo., 3mo., 6mo. and lyr. implied volatilities. 
Our results indicate that the modulation of autoregressive trend following us- 
ing derivative-based signals significantly improves the fit to the distribution of 
times between successive sign flips in the underlier time series. 



1. Introduction 

The spot market for major currencies is by far the most liquid financial market in 
the world. It is also notorious for its speculative nature, with over a trillion dollars 
changing hands daily. Most of this transaction volume is driven by endogenous 
'herd' dynamics, either technical or behavioral in nature. Over the past decade, the 
progressive abundance of publically available data at increasingly higher frequencies 
have made the FX market a testing ground for innovative models that probe where 
general equilibrium economics left off, beyond the Efficient Markets Hypothesis. 
On the one hand, empirically derived 'stylized facts' have become widely quoted 
as hallmarks of this new frontier. These include statistically challenging estimates 
of nonlinear, dynamics behavior, like clustered volatility (i.e. significantly positive 
serially correlation of realized volatility), long memory (i.e. subexponential decay 
of return autocorrelation function) and fat tails (i.e. probablity of extreme returns 
far exceeding that implied by the normal) PIE]. 

On the other hand, modelers have been moving their focus from models based 
on homogeneous, autonomous, sophisticated optimizing rational agents to hetero- 
geneous, strongly coupled agents with explicitly limited rationality. This represents 
a fundamental paradigm shift, necessitating new analytical tools to answer innova- 
tive questions. One modeling direction involves agents that choose among a small 
set of standard strategies. The most commonly studied trade-off is between trend 
following and contrarian trading. In this context, the question is how to optimally 
time the agents' strategic switching [Tl 131 ITT1 . 

Various attempts have been made to modulate pure trend following jSJ . Our 
approach extends earlier work in using multiplicative signals 0] to adjust endoge- 
nously the degree and timing of trend following. In this paper we present a hybrid 
dynamical system in which the continuous trend following process is multiplied by 
a volatility-driven signal |S] . We exploit derivative market information to construct 



Date: May 15, 2006. 



1 



2 



TED THEODOSOPOULOS AND ALEX TRIFUNOVIC 



Boolean logic rules that instantiate the desired switching between trend-following 
and contrarian strategies |%1 113). 

Our process consists of an AR(n) model multiplied by the symbolic dynamics 
of the term structure of at-the-money (ATM) implied volatility, which is assumed 
to diffuse lognormally. This simple model replicates the fat tails of the return 
distribution, the positively autocorrelated realized volatility and and the geometric 
phase plot pattern which characterise the actual data but are often absent in pure 
trend following simulations. Moreover, we are able to show a particularly accurate 
fit to the right tail of the 'inter-slip' times, i.e. the times between successive strategic 
switches from a follower to a contrarian rule. 

2. The Hybrid Model 

Let X k represent the price return time series, Y^ k the ATM volatility for the 
j th term at time k and e k an iid standard normal sequence. Let h : (—00,00) x 
[0,oo) — > { — 1,0,1} be the Boolean function that encodes the desired symbolic 
dynamics. Here we will use h{x,y) — X(y,oo)( x ) ~ X(-oo,-y)( x )- Consider the three 
parameter family of hybrid systems given by: 



(2.1) X k = n l g k 2J X k -i + ae k 

z=l 

for k > n, where a is a control parameter and 

(m 
2Mii,*-i*k-iiO),i 

We will model the volatility process as a driftless m-dimensional correlated random 

2 

walk, i.e. there exists a positive definite, symmetric matrix C S lZ m such that for 
all t = (ti, t 2 , . . . , t m )' and for all k > 0, 



E 

and for k\ ^ k% 



exp { \/ T T^Z ( Y i- k ~ Y i,k-i)tj 



E [(Y., fel - Y. ifcl _ x ) (Y, fc2 - Y. M -i)] = 0. 

We have simulated 10 6 steps of the hybrid model using the historically deter- 
mined covariance matrix for the ATM volatility term structure C We have used 
m = 4 for all simulations in this paper. Figure^shows the comparison of the hybrid 
model (12.11) on the far right column, against the pure trend following model 1 in the 
far left column and the actual daily returns of the JPY/USD exchange rate for five 
years in the middle column (2171 data points). The first row shows the comparison 
of the normal probability plots. This comparison shows that pure trend following 
generates thin tails, while the hybrid model matches well the fat tails characterizing 
the actual return distribution. In both cases n = 2 was used to generate the graph. 

The bottom row shows the comparison of phase plots, i.e. the relationship 
between X k and X k ^\. Naturally the pure trend following model exhibits very 
highly serially correlated returns, particularly because we've kept the noise level 
very low (a = 10 -4 ). Observe that the actual phase plot deviates noticeably from 



^This is equivalent to 12.11 with = 1. 



HYBRID DYNAMICS FOR CURRENCY MODELING 3 



Pure Trend (n=2) Actual Daily Returns (JPY/USD) Hybrid Model (n=2) 




-0.05 0.05 -0.06-0.04-0.02 0.02 -1 1 




-0.1 0.1 -0.04 -0.02 0.02 0.04 -1 1 

X k-1 X k-1 X k-1 x10~ 3 

Figure 1. The hybrid model matches better the distributional 
properties of the actual time series than pure trend following. 



the circular pattern expected from a random walk. The hybrid model generates a 
more square-like phase plot pattern, that nevertheless matches the phase plot for 
the actual data much better than the linear pattern of the pure trend following 
model. 

We have also investigated the degree to which the hybrid model generates re- 
alized volatility patterns consistent with the clustering phenomenon document in 
numerous empirical studies. Figure [21 shows the three way comparison of the bi- 
weekly realized volatility phase plot. In all cases there is a strong positive correla- 
tion, which is consistent with the previously reported volatility clustering. Indeed 
the resulting pattern from the hybrid model simulation better matches the low 
dispersion observed in the actual data. 



3. Symbolic Dynamics 

We proceed now to investigate in more detail the effect of the symbolic dynamics. 
Clearly, the multiplicative modulation of the AR(n) model with gk serves to flip 
the sign of the simulated returns at random intervals, determined by the direction 
of movement of the implied volatility. 
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Figure 2. Volatility clustering can be seen from the robust posi- 
tive slope linear patterns in the phase plots of realized volatility. 



Let p Q .™i m denote the path measure induced by l|2. 1[) with parameters a, n and 
m. Of course 

m 

^2h(Yj^ k -Yj-, fc _i,0) e {-m,-m + 2,-m + 4, ...,m- 4,m- 2,m} 

3=1 

and therefore 



(3.1) 



q x = lim P oAm (X fc+1 X fc < 0) = Pr (g k+1 < -1) 



which is independent of k since the volatility process has independent increments. 

Similarly for n — 2 we have: 

(3.2) 

q 2 = lim P Q < 2 <™ ({X k+1 (X k + X fc _i) < 0} n {X fc X fe _a > 0}) = Pr ( 5fe+1 < -1) 

a— >0 

Thus, in general, 
(3.3) 



q n = lim P 

a-+0 



' l ( j-X'fcS^*-* < °| n I n { X k-iXk-i-i > 0}|J = Pr [g k < -1) 

Thus qi = q% = . . . = q n = q. The driftless nature of the volatility diffusion 
implies that the resulting measure on {— to, —to + 2, — m + 4, . . . , to — 4, to — 2, to} 
is symmetric. Thus, only [^1 unknowns need to be determined. Using l3~3f we see 
that 



(3.4) 



1 - Pr -1 < YJLi h (Y jt k ~ Y jtk - U 0) < 1 
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FIGURE 3. The hybrid model matches better the distributional 
properties of the actual time series than pure trend following. 



Let T = and 

Ti — min < k > Tg_i 



x k J2^k- t <o n f| {Jffc-iJCfc.i.i > 0} 



for £ > and Tg = T{Tt_\. This sequence of stopping times represents the periods 
while the process X k maintains a fixed sign. The independent increments property 
discussed above is sufficient in the case n — 1 to obtain 



(3.5) 



Pr (n > a) = (1 - q y 



Figure |21 shows the complementary cumulative density function of r. We observe 
that the slope of this graph provides a clear indication of the advantage of the 
hybrid dynamics presented here in modeling the empirical data compared to pure 
trend following. Specifically the simulation of the hybrid model with n = 2 gives an 
estimated exponential decay rate of —1.6 which is within 3% of the value obtained 
from the JPY/USD data. On the other hand, the pure trend following model leads 
to a much shallower exponential decay rate of —1. 

Using (|3.5() we would expect an exponential decay rate equal to log(l — q) for 
a hybrid model with n = 1. The exponential decay rate for a simulation of the 
hybrid model with n = 1 is around —0.6371 (data not shown) which would lead to 
q = 0.4712. The resulting measure on {—to, —to + 2, —to + 4, . . . , to — 4, to — 2, to} 
is highly bimodal, as would be expected due to the high positive correlation among 
the implied volatilities. 

When n > 1 (13.51) no longer holds because of the longer history involved in the 
rolling average. Instead of the geometric model which gave rise to (|3.5|) we can 
use a higher order negative binomial process to accommodated the higher values of 



6 



TED THEODOSOPOULOS AND ALEX TRIFUNOVIC 



n. In particular, for a general n, we will approximate the distribution of re with a 
negative binomial with parameters q and n: 



The simulation which generated the data on the right-most panel in figure used 
n = 2. In this case, using l|3.6[l we would expect an exponential decay rate approx- 
imately equal to log(l — q) — 1, where we have used logo; = x for small x. Thus, 
the observed exponential decay of —1.6 corresponds to q = 0.4512. This value is 
slightly lower than that for the case n = 1 as expected because the averaging of two 
consecutive terms makes it somewhat more difficult to flip signs (i.e. the probability 
of long periods with the same sign decays slightly more slowly). 



We have described a simple hybrid dynamical system that models the modulation 
of trend following strategies for exchange rates using information from the market 
of currency derivatives. The resulting family of models was shown to better capture 
distributional characteristics of the empirical return distribution for the JPY/USD 
exchange rate. 

The particular form of the symbolic dynamics we used can be interpreted as 
a qualified majority of the implied volatility signals. Clearly, different Boolean 
expressions would lead to different patterns in the distribution of random times Tg 
at which the process changes sign. It is desirable to classify symbolic dynamics 
with respect to their effect on the resulting distribution of 'inter-flip' times. 

So far we have used a driftless model of the volatility term structure. Adding a 
true second moment signal (e.g. volatility drift) and an explicit third moment signal 
(e.g. risk reversal) can serve to make the symbolic modulation of trend following 
more nuanced. It would be useful to develop analytical tools, like the exponential 
decay rate approximation described in section [3J to detect the effect of different 
derivative signals in the empirically observable distributional patterns of inter-flip 
times. 
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